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ABSTRACT 

Characterization of microlensing planets requires modeling of observed light 
curves including many parameters. Studying the dependency of the pattern of 
light curves on the lensing parameters and the correlations between the parame- 
ters is important to understand how the uncertainties of the planetary parameters 
are propagated from other parameters. In this paper, we show that despite the 
apparent complexity of the pattern of light curves of planetary lensing events, 
the correlations between the lensing parameters can be understood by studying 
how the parameters affect the characteristics of lensing light curves such as the 
height and width, the caustic-crossing time scale, and the location and duration 
of planetary perturbations. Based on analytic arguments about the dependency 
of light curve features on the parameters, we obtain the correlations for the two 
representative cases of planetary events. We also demonstrate the applicability of 
the correlations to general planetary events by actually obtaining the correlations 
from modelings of light curves produced by simulations. 

Subject headings: gravitational lensing 



1. Introduction 



A microlensing event occurs when an astronomical object (lens) is closely aligned with 
the line of sight toward a background star (source). Microlensing causes change of the 
source star bright ness and the resulting light curve is characterized by its smooth variation 
( iPaczyhski I Il986l ). If the lensing object is a star and it contains a planet, the resulting 
light curve can exhibit a discontinuous signature of the planet on the smooth light curve 
of the primary-ind uced event, and thus mic r olensing can be used as a method to search for 
extrasolar planets ( IMao. fc Paczynskilll99ll ; iGould &: Loebl Il992l ) . Microlensing is sensitive 
to planets that are generally inaccessible to other methods, in particular cool planets at or 
beyond the snow line, very low-mass planets, planets orbiting low-mass stars, free-floating 



-2- 



planets, and even planets in external galaxies. Therefore, when combined with the results 
from other surveys, microlensing planet searches can yield an accurate and complete census 
of the frequency and properties of planets. Since the first discovery in 2004, 9 microlens- 



ing planets 


lave been reported f 
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2010|; 
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2010). 



Characterization of microlensing planets requires modeling of observed light curves. 
This modeling process requires to include many parameters because the pattern of lensing 
light curves and the signals of planets take different forms depending on the combination 
of these parameters. Therefore, studying the dependency of the pattern of light curves on 
the lensing parameters and the correlations between the parameters help to understand how 
the uncertainties of the planetary parameters are propagated from other lensing parameters. 
This also helps to establish observational strategies for better characterization of planets. 
However, it appears that the correlations between the lensing parameters are very complex 
due to the enormous diversity of lensing light curves resulting from the combinations of the 
numerous parameters. 

In this paper, we show that despite the apparent complexity of the pattern of light 
curves of planetary lensing events, the correlations between the lensing parameters can be 
understood based on the dependency of the characteristic features of lensing light curves on 
the parameters. We provide the correlations for the two representative cases of planetary 
events. We also demonstrate the applicability of the correlations to general planetary lensing 
events by actually obtaining the correlations from modelings of light curves produced by 
simulations. 



2. Parameters 

The microlensing signal of a planet is a brief perturbation to the smooth standard 
light curve of the primary-induced single- lensing event. Therefore, the parameters needed to 
describe planetary lensing light curves are broadly divided into two categories. The first set 
of parameters is needed to describe the light curve of a standard single-lens event produced 
by the star hosting the planet. These parameters include the closest lens-source separation 
normalized by the Einstein radius, (3 (impact parameter), the time of the closest lens-source 
approach, to, the time required for the source to transit the Einstein radius of the lens, £e 
(Einstein time scale), the flux from the lensed star, F$, and the blended flux, Fb- With these 
parameters, the single-lensing light curve is represented by 

F(t) = F s A(t) + F B , (1) 
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Here w represents the lens-source separation normalized by the Einstein radius 9 E . These 
parameters characterize the global shape of lensing light curves such as the height and width. 

Besides the single-lensing parameters, additional parameters are needed to described the 
detailed structure of the perturbations induced by planets. These parameters include the 
planet/star mass ratio, q, the projected star-planet separation normalized by the Einstein 
radius, s, and the angle between the source trajectory and the star-planet axis, a. Since 
planetary perturbations are produced in most cases by close approaches or crosses of source 
stars over caustics, an additional parameter of the source size normalized by the Einstein 
radius, p*, is needed to describe the deviation of the perturbation affected by the finite-source 
effect. From the combinations of the lensing parameters, the light curves of planetary events 
exhibit various patterns. 



3. Correlations 

3.1. Single-Lensing Parameters 

The correlations between the single-lensing parameters can be found by investigating 
how the global features of lensing light curves such as the height and width vary depending 
on the parameters. The height of the light curve is determined by the combination of the 
impact parameter (3 and the blended light ratio, /b = Fb/F$. When affected by blended 
flux, the impact parameter estimated based on the apparent height of the light curve is 
related to the blended light ratio as 

ft = [2(1 - A m 2 ax , B r 1/2 - 2] 1/2 ; A max , B = (3) 

where A maxB is the apparent peak magnification. Then, with the increase of /b, A maxB 
decreases, and thus (3b increases, implying that as the blended light increases, the peak mag- 
nification appears to be lower, and the resulting impact parameter becomes larger. Therefore, 
it is found that the blended flux ratio /b and the impact parameter (3 are correlated. 

Blending affects not only the height but also the width of light curves. The event time 
scale estimated from the width of the blended light curve, £e,b, differs from the true value, 
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t&. The relation between the two time scales is (jHanlll999l ) 



iE ^H1t^J • (4) 

Here /3 t h,B and /3 t h represent the threshold impact parameters of the source trajectory for 
event detections with the presence and absence of blended flux, respectively. Because of 
blending, the threshold magnification is increased by Ah,B = Ah(l + /b) — /b and thus the 
corresponding threshold impact parameter is lowered by /3 t h,B = [2(1 — A^b) 2 — 2] 1//2 , 
implying that to be detected the source star of a blended event should approach the lens 
closer than an unblended event. From equation (j3J, it is found that the increase of blended 
flux causes the measured event time scale shorter. Therefore, the blended flux ratio /b and 
the Einstein time scale t^ are anti- correlated. Since the blended flux ratio is correlated with 
the impact parameter, while it is anti-correlated with the time scale, it is found that the 
impact parameter (3 and the time scale t E are anti-correlated. 



3.2. Binary Parameters 



The major constraint on the the normalized planetary separation s comes from the 
location of the planetary perturbation on the light curve. This is because the location of the 
caustic with respect to t he primary lens, u p , is re lated to the position vector to the planet 
from the primary, s, by (IGriest fc Safizadehlll998l ) 

1 



(5) 



The primary-caustic separation, u p , varies differently depending on whether the planet is 
inside (s < 1) or outside (s > 1) of the Einstein radius. For planets with s > 1, the 
caustic separation u p increases with the increase of the planetary separation s. By contrast, 
the caustic separation decreases as the planet separation increases. In addition, the type 
of the perturbation differs depending on the location of the planet. When the planet is 
located outside the Einstein ri ng (s > 1), it perturbs the major image among the two images 
produced by the primary lens (IGaudi fc Gouldlll997l ). resulting in higher magnifications than 
those of the adjacent part of the single- lensing light curve during the perturbation. When 
the planet is inside the Einstein ring (s < 1), by contrast, it perturbs the minor image and 
the magnifications during the perturbation are lower than those of the single-lensing light 
curve. From the relation in equation ([5]) combined with the geometry of the lens system 
presented in Figure (TJ it is found that the relative location of the perturbation on the light 
curve is related to the time of perturbations, t p , and other lensing parameters by 

t P — to 



tE 



uc 



2x1/2 



(6) 



-5 - 



Then, for planets with s > 1, as s increases, u p increases, and thus the impact parameter 
(3 should increase to match the constraint given by the location of the perturbation on the 
observed light curve. For planets with s < 1, by contrast, u p decreases with the increase of s, 
and thus (3 decreases. Therefore, it is found that the planetary separation s and the impact 
parameter (3 are correlated for planetary events with major-image perturbations (s > 1) and 
anti- correlated for events with minor-image perturbations (s < 1). 

The correlation between the planetary separation s and the mass ratio q is mainly given 
by the duration of the planetary perturbation. The duration is proportional to the size of 
the caustic. The size of the cau stic as measured by the width along the primary-planet axis 



is related to s and q ( IHanl 120061 ) by 



f 4^/^(1 _ s -2)] for s > x, 
c \2q 1 / 2 (K-K- 1 + KS - 2 ) fors<l, K) 

where k 2 = [cos 29 ± (s 4 - sin 2 29f/ 2 }/{s 2 - s~ 2 ) and 9 = vr/2 ± sm _1 (V3s 2 /2)/2. In the 
limiting cases of planetary separations of s 3> 1 and s<l, the relations are expressed in 
compact forms of 

j 4g i/2 s -2 for s » 1, 

Wc ~ UVZq 1/2 s*/A for S «l, { ) 

respectively. Then, from the relation between s and q combined with the constraint of the 
perturbation duration measured from the light curve, it is found that the planetary separation 
and the mass ratio should be correlated for planetary events with major-images perturbations 
(s > 1 ) and anti- correlated for events with minor-image perturbations (s < 1). 

The main constraint on the normalized source radius p* is provided by the duration of 
caustic-crossings (caustic-crossing time scale). The caustic-crossing time scale is related to 
the normalized source radius p+ and the Einstein time scale t& by 

Ai cc = -A-t E , (9) 
smip 

where ip represents the angle between the caustic and the source trajectory. Then, to match 
the constraint of the caustic-crossing time scale measured from the observed light curve, the 
normalized source radius p* and the Einstein time scale should be anti- correlated. 

In Table (TJ we summarize the correlations between the lensing parameters. Here, the 
correlations between the pairs of parameters not mentioned in the text are deduced from 
the correlations with other parameters. For example, we deduce the correlation between the 
normalized source radius p* and the impact parameter (3 based on the correlations between 
P*-£e an d t-£-f3 parameter pairs. In the table, we mark for the pairs of parameters that 
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are correlated, while the correlation is marked by "\" for the pairs of parameters that are 
anti-correlated. 



Gould fc Welch! ( 119961 ) studied the relation between the normaliz ed source radius and the 



impact parameter for high-magnification single-lens events. See also iGaudi fc Gouldl (119971 ) 
for the relation between the lens-source proper motion, //, and the time scale, although fi is 
not a standard planetary lensing parameter. 



4. Demonstration 

In the previous section, we investigated the correlations between the lensing parameters 
based on analytic arguments about the dependency of the characteristic features of lensing 
light curves on the parameters. In this section, we demonstrate that the correlations are 
applicable to general planetary microlensing events by actually obtaining the correlations 
from modelings of light curves produced by simulations. 

We produce two light curves of planetary lensing events, where the individual curves 
represent those of events with major and minor-image perturbations, respectively. The light 
curves are produced considering the strategy of the current planetary lensing experiments 
where events are detected through modest-cadence survey observations and perturbations 
are densely covered by follow-up observations. We set the photometric uncertainty by as- 
suming that the photometry follows photon statistics with a 1% systematic uncertainty and 
the deviations of the data points are Gaussian distributed. The other factors affecting the 
photometry such as the source brightness and blending are based on the values of typical 
Galactic bulge events. Figure [2] shows the light curves of the events produced by the simula- 
tion where the upper and lower panels are those of the events with major and minor-image 
perturbations, respectively. The inset in each panel shows the geometry of the event where 
the straight line with an arrow is the source trajectory and the temperature scale represents 
magnifications where a brighter tone implies a higher magnification. 

We search for the solution of the lensing parameters by conducting modeling of the 
simulated light curves. In modeling planetary lensing light curves, it is difficult to conduct 
brute-force search throughout all parameter space due to the large number of parameters. 
In addition, it is difficult to conduct simple downhill approach due to the complexity of the 
X 2 surface. We, therefore, use a hybrid approach where grid searches are conducted over 
the space of parameters of s, q, and a and the remaining parameters are allowed to vary 
so that the model light curve results in minimum \ 2 at each grid point. We use a Markov 
Chain Monte Carlo method for \ 2 minimization. Once the solutions of the individual grid 
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points are determined, the best-fit model is obtained by comparing the x 2 minima of the 
individual grid points. The uncertainties of the best-fit parameters are estimated based on 
the chains of solutions produced by modeling. For the computations of lensing magnifications 
including finite-source effect, we use the ray-shooting method where a large number of rays 
are shot from the image plane, bent according to the lens equation, and arrive on the 
source plane. Then, the magnification corresponding to the location of a source star is 
computed by comparing the number density of rays on the source star with the density on 
the image plane. We minimize the computation time by restricting the region of ray shooting 
around the images and thus minimizing the number of rays needed for the computations 
of finite-source magnifications. We als o apply semi-analytic hexadecapole approximation 



(IPejcha fc Heyrovskyil2009l ; lGouldll2008l ) for magnification computations in the region where 



the effect of source size is not important. 

In Figure [3] and HI we present the results of the modeling for the planetary events with 
major and minor image perturbations, respectively. The individual panels of each figure 
show the contour plots of x 2 m t ne spaces of the combinations of the lensing parameters. To 
directly compare the correlations, we arrange the panels according to the same order of the 
parameters presented in Table [TJ From the comparison, it is found that the correlations found 
from modeling coincide with those predicted based on analytic arguments. This implies that 
the correlations presented in Tabled] are applicable to general planetary microlensing events. 



5. Conclusion and Discussion 

We investigated the correlations between the parameters of planetary lensing events. 
From this, we found that the correlations could be understood by studying how the lensing 
parameters affect the characteristics of lensing light curves such as the height and width of the 
light curve, the caustic-crossing time scale, and the location and duration of perturbations. 
Based on analytic arguments about the dependency of the features of lensing light curves on 
the parameters, we obtain the correlations. We also demonstrated the applicability of the 
correlations to general planetary lensing events by actually obtaining the correlations from 
modelings of light curves produced by simulations. 

Understanding the correlations between lensing parameters can help to setup obser- 
vational strategies for better constraints of the planetary parameters. For example, the 
correlations between the blending and planetary parameters imply that blending affects de- 
terminations of the planetary parameters and thus de-blending or precise determination of 
the blending parameter is important to better constrain the planetary parameters. Several 
methods can be applied to resolve the blending problem. Photometrically, it is known that 
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good coverage of the peak region an d wings of light curves by follow-up observations help 
to constrain the ble nding para meter (jThomas fc Griestll2006l ). High resolution ima ging; from 
space observations (|Hanlll997l ) or ground-based AO imaging (IBennett et al.ll2010f ) can help 
to resolve source stars from blended stars. Astrometric measurement of the blended imag e 
centroid can also help to identify the lensed source among blended stars (lAlard et al.lll995l ). 
Precise and dense coverage of the perturbation is another way to constrain the planet param- 
eters. This is because not only the location of perturbations but also their shape provides 
constraints on planetary parameters^] Then, even if the location of the perturbation is un- 
certain due to severe blending, it is still possible to constrain the planetary parameters from 
the shape of the perturbation. 



This work is supported by Creative Research Initiative program (2009-0081561) of Na- 
tional research Foundation of Korea. 
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Fig. 1. — The geometry of the lens system showing the relation between the location of 
the planetary perturbation and lensing parameters. The large dashed circle represents the 
Einstein ring and the small closed figure connected by concave curves is the caustic induced 
by the planet. 
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Table 1. Correlations Between Lensing Parameters 
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Fig. 2. — Light curves of planetary lensing events with major (upper panel) and minor- 
image (lower panel) perturbations. The inset in each panel shows the geometry of the 
event where the locations of the primary and planet, source trajectory (straight line with an 
arrow), caustic (small closed figure), and Einstein ring (big dashed circle) are marked. The 
correlations between the lensing parameters of the individual events are presented in Fig. [3] 
andHl respectively. 
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Fig. 3. — Contours of x 2 i n the spaces of various combinations of lensing parameters for 
the planetary event with a major-image perturbations. The result is based on the modeling 
of the light curve presented in the upper panel of Fig. [2j For direct comparison of the 
correlations to those presented in Table [U, the order of parameters are arranged accordingly. 



-14- 




Fig. 4. — Same as in Fig. [3] except that the correlations are for the planetary event with a 
minor-image perturbation. 



